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Coastal  subsidence  increases  flooding  in  low  lying  coastal  regions. 
Moreover,  it  disturbs  the  equilibrium  profile,  and  allows  waves  to  erode 
bluffs  formerly  above  the  reach  of  wave  uprush.  Ensuing  adjustment  of  the 
profile  drives  the  shoreline  farther  landward.  Guidance  is  needed  for 
obtaining  quantitative  estimates  of  the  shore's  response. 

The  mean  surface  elevation  of  Lake  Michigan  rose  0.5  meter  during  a recent 
4-year  period.  Concurrently,  major  elements  of  the  submerged  profile  — ^ 
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responded  by  building  upward  and  migrating  26  meters  landward.  Approximately  8 
meters  of  beach  were  lost  due  to  submergence  beneath  the  elevated  lake  surface; 
and  an  additional  6 to  7 meters  were  lost  due  to  erosion.  The  shoreline, 
however,  lagged  beliind  the  rest  of  the  profile  in  adjusting  to  the  higher  water 
levels . ^ 


Over  a lon^er^  period,  certain  sections  of  Lake  Michigan  have  undergone 
relative  subsidence  as  a consequence  of  broad  regional  tilting  of  the  earth's 
crust.  It  is  estimated  tliat  during  the  last  century  the  Lake  Michigan  basin 
tilted  0.06  to  0.09  meter  per  100  kilometers  along  its  axis.  Shore  recession 
over  the  last  century  increased  at  a rate  of  19  ± 10  meters  per  100  kilometers 
in  the  direction  of  greater  subsidence. 

Other  coastal  areas  with  similar  geomorphology  and  wave  exposure  can  be 
expected  to  recede  at  rates  similar  to  those  indicated  above  if  subjected 
to  the  same  subsidence.  The  initial  response  to  rapid  subsidence  may  be  on 
the  order  of  50  units  of  retreat  for  each  unit  of  subsidence.  Profile  retreat 
is,  however,  a nonlinear,  time-dependent  function  of  subsidence,  and  for  slower 
subsidence,  shore-eroded  sediments  become  spread  over  a broader  profile,  pro- 
ducing a larger  ratio  of  shore  retreat  per  unit  of  subsidence.  Estimates  of 
shore  recession  due  to  slow  crustal  motion  of  the  Great  Lakes  basin  indicate 
response  ratios  between  120:1  and  390:1.  Furthermore,  according  to  the  concept 
of  mass  balance,  the  long-term  response  ratio  should  also  depend  on  the  volume 
and  size  distribution  of  sediments  being  supplied  to  the  nearshore  profile  by 
shore  erosion.  The  lakeshore  response  was  found  to  increase  several  fold 
where  inadequate  backshore  deposits  supplied  less  beach  materia-  per  unit 
of  recession. 
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Cooi  tai  iub-i  idenccL  incfie-cau  ^coding  -tit  tow  lying  cx>aitat  ft£.gioni, 
McAcouc-i,  a.  dtitiLibi  the  eguitibiium  pA/j^ile,  and  attouii  wavei  to  etude 
6^u(5i5i  {ot/netlij  above  the  teach  of  wave  uptush,  Emuing  adjuAtment  the 
pte  ^ite  dtevei  the  ihoteti-iie  ^a-ithet  taiideoAd.  Guidance  it  needed  fot 
obtainiitg  quantitative  eAtimatei,  o(^  the  Ahote’A  teApoiiAe. 

The  mean  AuA^ace  etevatcon  Labe  Michigan  toAe  0,5  m dating  a tecent 
(puA  yeat  petiod.  ConcutAcnUy,  mx/ot  elementA  05  tlic  Aubmetged  ptaiile 
tciponded  by  building  upevatd  and  migtating  16  m landwatd.  Approximately  t 
m oi5  beaeh  wete  loAt  due  to  Aubmetgence  beneath  the  elevated  lake  Aut^ace; 
and  an  additional  6 to  7 m wete  loAt  due  to  The  Ahoteline,  howevet, 

tagged  behind  the  teAt  oi  the  profile  in  adjuAting  to  the  highet  loatet 
tevelA . 

Over  a longer  period,  certain  AecteonA  0^  Lake  Michigojx  tuxve  u}ider~ 
gone  relative  AubAidence  oa  a conAequenee  o^  broad  regional  tilting  oi  the 
earth 'a  cruA  t.  It  La  eAtimated  that  during  tiie  laAt  century  the  Lake 
Michigan  baAin  tilted  0.06  to  0.09  m per  WO  kilometetA  along  itA  axlA. 

Shore  reccA  A ion  over  the  laAt  century  increaAed  at  a rate  19  i.  W m per 
100  kilometerA  in  the  direction  oi  greater  AubAidence. 

Other  cocLAtal  areoA  with  Aimilar  geomorphology  and  luave  expcAure  can 
be  expected  to  recede  at  ratcA  Aimilar  to  thoAe  indicated  above  ii  Aubjected 
to  tlie  6ame  iubii-dence.  The  iytitiaZ  rcApoiiAe  to  rapid  AubAidence  may  be  on 
the  order  oi  50  unitA  oi  -retreat  ior  each  unit  oi  AubAidence.  Proiile  re- 
treat iA,  however,  a non-linear,  time -depen  deni  (unction  oi  AubAidence,  and 
ior  Alower  AubAidence,  Ahore-eroded  AcdijnentA  become  Apread  over  a broader 
pro^e,  producing  a larger  fjotio  oi  Ahore  retreat  per  unit  oi  AubAidence. 
EAtiinatcA  oi  Ahore  reccAAion  due  to  Alow  eruAtal  motion  oi  the  Great  Lake  a 
baAin  indicate  reApouAc  ratioA  between  120:1  and  390:  1.  furthermore, 
according  to  the  concept  oi  moAA  balance,  the  Cong  term  reAponAc  ratio 
Ahould  oLao  depend  on  the  volume  aird  Ai  .e  diAtribution  oi  AedimentA  be-in-^ 
Aupplled  to  the  nearshore  proiile  by  shore  en^sion.  The  lakeshore  responAe 
ivas  ioiuid  to  increase  several  iold  where  inadequate  backshore  depoAitA 
supplied  less  beach  material  per  unit  oi  recession. 

I INTRODUCTION 

1 . Aim  and  Approach 

Failure  to  consider  the  effects  of  subsidence  may  result  in  a serious 
error  in  projection  of  long  term  shore  retreat  at  certain  localities.  Even 
moderate  rates  of  subsidence  can  increase  coastal  erosion  and  submerge  low 
lying  coastal  regions.  Furthermore  many  conditions  brought  about  by 
development  of  shore  property,  e.g.,  lowering  of  the  water  table  accompany- 
ing ground  water  withdrawal,  soil  compaction  in  response  to  Increased  sur- 
face loads,  vibrations  associated  with  construction  operations  and  later 
traffic,  can,  as  documented  by  other  papers  in  this  symposium,  contribute  to 
subsidence;  and  thereby  accelerate  shore  retreat.  E.stimates  of  potential 
shore  erosion  should,  therefore,  include  allowances  for  the  continuation  of 
documented  subsidence,  and  even  for  ejcpected  Increases  in  the  subsidence 
rate  when  the  state-of-knowledge  permits  such  predictions,  based  on  develop- 
ment factors.  DISTRIBUTION  STATEMliNT  A i 
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No  guidance  is  available  at  present  however,  for  estimating  the  amount 
of  erosion  which  will  result  from  a given  rate  of  subsidence.  More  funda- 
mental questions  about  coastal  currents,  waves,  and  the  interactions  of 
fluid  and  sediment,  have  yet  to  yield  solutions;  and  In  the  absence  of  a 
theoretical  understanding,  empirical  correlations  between  shore  retreat  and 
causative  factors  have  assumed  subsidence  per  se  can  be  Ignored  (Caldwell, 
1959;  Richardson,  1976;  and  others).  Usually  the  effects  of  subsidence  on 
shore  retreat  rates  are  relatively  small  compared  with  the  total  variation 
in  retreat  rates.  Moreover  even  where  subsidence  has  a significant  cumula- 
tive effect  the  shore  damage  will  tend  to  be  felt  during  the  storm  events. 
Thus  the  present  uncertainty  concerning  the  relationship  between  coastal 
subsidence  and  shore  retreat  should  not  be  surprising. 

Because  subsidence  js  a minor  variable  among  many  affecting  the  rates 
of  shore  erosion,  resolution  of  its  contribution  will  depend  on  shore 
changes  measured  at  sites  of  extreme  subsidence,  or  over  long  periods  of 
time.  The  purpose  of  this  paper  is  to  estimate  the  effects  of  coastal  sub- 
sidence using  first,  data  on  Lake  Michigan  shore  retreat  during  four  years 
of  rapidly  rising  take  levels;  and  secondly,  using  historic  data  on  the  120 
year  retreat  rate  along  sections  of  the  Lake  experiencing  different  rates 
of  relative  subsidence. 


2 . Terminology. 

For  the  purpose  of  this  paper,  precise  definitions  have  been  formu- 
lated, that  refine  the  meaning  of  several  familiar  words.  SufameAgencc  will 
refer  to  the  sinking  of  a coastal  area  relative  to  the  mean  water  surface. 
Submergence  may  result  from  either  land  subsidence  or  elevation  of  the  wate: 
surface.  Emergence  refers  to  the  opposite  relative  displacement,  and  when 
expressed  numerically,  both  will  refer  to  lengths  (L)  measured  in  the  verti- 
cal. Coastal  planners  and  property  owners  are  often  more  interested  in  the 
resulting  horizontal  change  in  shoreline  position.  Tfiaiiig'izi6^Lon  (ant. 
•'legTte-li-c.on)  will  refer  to  the  horizontal  distance  that  the  shoreline  moves 
in  direct  response  to  submergence  (emergence).  The  ilWittA-iie.  is  the  inter- 
section of  the  beach  with  the  mean  water  surface,  or  some  other  specified 
datum  such  as  low  water  (LWD)  or  mean  high  water  (MHW) . The  shoreline 
divides  the  beach  into  iho>ie.  (subaerial)  and  nzd-uhofid.  (submerged)  zones. 

Total  lateral  migration  of  the  shoreline  can  be  more  or  less  than 
transgression  (regression)  depending  on  whether  erosion  or  deposition  pre- 
vails at  the  shoreline.  Deposition  refers  to  the  accumulation  of  material 

(ilO 


on  a surface  (M/L^  or  L) ; erosion  refers  to  its  removal.  The  lateral 
migration  of  a specific  contour  will  be  referred  to  as  p^Q>uidation  (L)  if 
the  contour  moves  toward  the  center  of  the  basin  and  as  fizcuiion  if  the 
contour  moves  away  from  the  basin.  Sho'ietcm  fieXKZOt  (Fig.  2)  will  be  the 
inclusive  term  referring  to  the  total  landward  horizontal  shift  or  the 
algebraic  sum  of  tAanig^ei6'ion  (a  function  of  submergence)  and  /lecfcAAcon 
(a  function  of  subsequent  erosion) . 


Figure  2.  Terminology  of  Retreat 
Transgression  = Az  cot  a 
Retreat  - transgression  + recession 

Shoreline  retreat  implies  that  there  has  been  either  local  recession 
or  transgression,  but  is  unspecific  as  to  which  (or  whether  both)  caused 
the  landward  shift  in  shoreline  position.  Figure  3 Illustrates  the  meaning 
and  heirarchy  of  the  discussed  terms. 
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Figure  3a.  Terminology  of  Figure  3b.  Terminology  of 

Vertical  Processes  Horizontal  Migration 

II  USE  OF  LAKE  LEVEL  FLUCTUATIONS  TO  MODEL  SUBSIDENCE 
1 . Lake  Level  Fluctuations. 

In  terms  of  submergence  and  erosion,  it  is  immaterial  whether  the  land 


Is  sinking  or  the  sea  rising.  Thus  a stable  or  even  uplifted  shore  can 
serve  as  a model  If  only  the  sea  Is  rising  facter  than  the  land.  Due  to 
their  distinctive  hydrologic  cycles,  the  Great  Lakes  can  be  used  to  model 
effects  of  subsidence  on  shoie  retreat. 

The  mean  surface  elevation  of  Lake  Michigan  has  a definite  annual 
cycle  averaging  about  IS  cm  In  amplitude.  This  seasonal  cycle  Is  super- 
imposed on  iong  teA/n  ^uctuaCioni  which,  though  not  so  regular  as  the 
seasonal  cycle,  are  nevertheless  clearly  evident  In  the  115-year  hydrograph 
of  annual  means  (Fig.  A) . 


Figure  4.  Annual  Mean  Surface  Elevations  of 
Lake  Michigan  (data  from  NOAA  records) 


Seasonal  and  long  term  lake  level  fluctuations  are  expressions  of 
meterologlc  and  climatic  variations.  As  shown  In  Figure  4,  the  annual  mean 
lake  level  generally  Increases  for  several  years  In  succession  before  going 
Into  periods  of  decline.  The  cumulative  effect  has  often  been  a consider- 
able shift  In  annual  mean  lake  levels,  e.g.,  more  than  a meter  between  1926 
and  1929;  0.83  meters  between  1950  and  1952;  and  1.45  meters  between  1964 
and  1973.  These  long  term  lake  level  fluctuations  are  much  more  rapid  (34, 
41,  and  16  cm/yr  respectively)  than  trends  In  sea  level  which  are  typically 
on  the  scale  several  mm/yr.  Runs  of  Increasing  lake  level  have  rates  more 
comparable  to  changes  In  relative  sea  level  at  sites  of  extreme  local  sub- 
sidence. This  suggests  using  measurements  made  to  determine  the  shore 
response  to  periods  of  rising  lake  levels,  as  Input  to  a model  for  esti- 
mating the  effects  of  coastal  subsidence. 

2 . Historic  Data  on  Bluff  Recession. 

Land  surveys  and  aerial  photographs  have  been  used  to  document  Lake 
Michigan  bluff  recession  In  excess  of  300  meters  (1116  feet)  In  121  years 
(Powers,  1958)  . Average  bluff  recession  rates  shown  In  figure  5,  represent 
nearly  1000  measurements  compiled  from  various  references  reporting  the 
position  of  bluffs  at  widely  scattered  sites  around  the  lake  shore.  The 
height  of  each  horizontal  bar  represents  the  recession  rate  obtained  by 
averaging  measurements  from  9 to  100  different  sites.  The  length  of  the  bar 
Indicates  the  various  time  spans,  from  2 to  127  years,  for  which  recession 
was  calculated.  In  this  format,  measurements  from  localities  having 
different  degrees  of  exposure  and  resistance  to  erosive  forces  are  lumped 
together,  but  two  Important  features  nevertheless  stand  out:  first,  the 
average  lakeshore  bluff  Is  evidently  subject  to  persistent  long  term  retreat. 


Secondly,  rates  of  retreat  are  not  uniform  through  time.  During  certain 
specific  short  periods,  average  recession  rates  increase  significantly 
throughout  the  lake. 

It  has  long  been  held  that  shore  erosion  was  most  rapid  during  periods 
of  above  average  lake  level.  Evidence  of  this  was  presented  by  Selbel 
(1972)  for  Lake  Michigan/Huron  and  by  Berg  and  Duane  (1968)  for  Lake  Erie. 
Some  investigators  nevertheless  still  express  doubt  as  to  the  validity  of 
this  relationship  (Larsen,  1972;  McKee,  1972)  . A comparison  of  the  rates 
of  recession  shown  in  Figure  5 with  the  mean  annual  lake  level  (shown  by 
the  continuous  jagged  line)  shows  direct  correlation  between  bluff  reces- 
sion rates  and  the  mean  elevation  of  the  lake  surface,  in  spite  of  the 
fact  that  the  data  collection  periods,  over  which  recession  rates  are 
averaged,  were  undesirably  long  and  not  ideally  suited  for  resolving  the 
effect  of  lake  level  changes.  Consider  for  example,  the  change  in  reces- 
sion at  IB  in  Figure  5.  This  average  rate  of  recession  is  based  on  AO 
point  measurements  along  a specific  length  of  shore.  A 500%  increase  in 
average  recession  rate  is  noted  between  1938-1950  and  1950-1955.  The  rate 
then  decreased  considerably  as  lake  levels  fell  in  the  late  50's.  Both 
lake  level  and  recession  decreased  still  further  in  the  early  sixties. 

Lake  levels  reached  a record  low  in  196A.  Recession  was  low  but  not  below 
the  long  term  average,  which  may  reflect  first  the  effect  of  available 
survey  intervals;  recession  and  lake  levels  were  still  moderately  high  in 
'60,  '61,  and  '62,  and  then  rose  again  rapidly  in  the  closing  years  of  that 
decade.  A second  possible  reason  for  above  average  recession  in  the  same 
decade  in  which  levels  dropped  to  a record  low,  is  that  the  studies  from 
which  the  IB  rates  (as  well  as  most  of  the  other  post  AO's  data)  were  taken, 
concentrated  on  the  more  critically  eroding  sections  of  the  lake. 

An  estimate  for  the  overall  long  term  recession  would  be  0.37  m/yr 
(1.2  ft/yr)  based  on  measurements  at  9A  stations  selected  from  data 
presented  by  Powers  (1958).  Powers'  data  provides  a good  estimate  of  the 
overall  historic  recession  because  Powers  chose  his  sites  systematically, 
and  the  9A  sites  used  here  were  all  originally  surveyed  between  1830  and 
1838.  Powers  determined  bluff  recession  by  surveying,  in  1956  and  1957, 
the  distance  of  the  shore  bluff  from  township  and  range  section-comers 
within  0.5  mile  of  the  lake.  Comparing  his  measurements  with  original 
government  surveys,  he  found  that  the  bluff  had  advanced  at  six  of  the  13A 
sites  (average  rate  of  advance  0.5  m/yr);  and  showed  no  change  at  four 
sites. 

3.  Profile  Adjustment  to  a Single  Period  of  Rapidly  Rising  Lake  Levels. 

The  response  of  the  beach  to  the  most  recently  rising  lake  levels  has 
been  monitored  at  six  stations  in  the  vicinity  of  Pentwater  Harbor  about 
midway  up  the  eastern  shore  of  Lake  Michigan.  Study  of  profile  changes  has 
provided  an  estimate  of  the  increase  in  shoreline  retreat  due  to  high  lake 
levels,  permitted  the  resolution  of  transgression  and  recession,  and  re- 
vealed changes  across  a broad  submerged  section  of  the  profile  which  are 
also  related  to  submergence.  The  dates  of  the  four  field  seasons  during 
which  the  six  stations  were  reprofiled,  together  with  the  change  in  lake 
level  between  the  field  seasons  (based  on  average  dally  means),  and  the 
mean  monthly  elevations  during  intervening  periods  are  given  in  Figure  6. 

a.  Shore . The  net  retreat  of  the  shoreline  over  the  four  year  study 
period  is  shown  in  Table  1.  In  spite  of  a slightly  higher  lake  level  in 
the  fall  of  1969,  the  shoreline  advanced  between  the  spring  and  fall  at  two 
of  the  six  profile  stations  (3  & 7)  because  a small  coastal  bar  merged  with 
the  shore.  Over  the  longer  period  from  spring  1969  to  1971,  a net  retreat 
developed  at  all  stations.  The  average  retreat  rate  for  the  two  year 
period  was  A m/yr,  but  there  was  still  a considerable,  random  variation, 
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in  retreat  among  the  different  stations.  Over  the  1967-1971  (45  month 
period,  longshore  variations  nearly  vanished  as  all  stations  approached  the 
average  retreat  rate  of  4 m/yr. 
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Figure  6.  Lake  Michigan  Hydrograph  Showing 
Differences  in  Mean  Lake  Levels  Between 
the  Various  Survey  Periods. 


Table  1.  Net  shoreline  retreat 


Station  number 

Spring 

to  l-all  1969 

Spring  1969  to 

1971  1967  to  1971 

3 

-1.5 

1.5 

13.4 

4 

1.5 

12.0 

16.7 

5 

2.5 

10.7 

14.6 

6 

5.3 

7.5 

15.2 

7 

-0.2 

5.8 

16.9 

8 

2.0 

12.5 

11.0 

Avg.  retreat  (m) 

1.3 

8.3 

14.6 

Coefficient  of 

variation  (m) 

1.4 

0.51 

0.  15 

Avg.  retreat 

rate  (m/yr) 

3.3 

4.1 

3.9 

A determination  of  the  exact  amount  of  >ie.CZii-Lon  depends  on  the  ele- 
vation where  the  measurement  is  made.  One  convenient  choice  is  at  the 
elevation  of  lake  surface  in  1967.  The  average  of  daily  means  during  the 
1967  field  season  was  176.30  m (International  Great  Lakes  Datum,  ICLD,  1961) 
and  the  average  for  the  whole  year  was  slightly  lower;  176.10.  The  positions 
of  these  two  ihofizlyinzi  as  well  as  the  LIllV  iho^eZ^e.  (175.80)  were  calculated 
No  datum  higher  than  the  176.1  could  be  selected,  however,  because  profiling 
was  terminated  at  water's  edge  in  1967.  If  a datum  lower  than  LWD  had  been 
selected,  it  would,  by  the  time  of  the  1971  survey,  have  Intersected  the 
profile  lakeward  of  a longshore  bar.  As  can  be  seen  in  Figure  7,  this  bar 
was  also  migrating  landward  as  lake  level  rose.  A coastal  bar  migrating 
landward  and  upward  can  cause  a sudden  anomalous  lakeward  advance  of  some 
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contours  on  an  otherwise  receding  shore.  Thus,  measured  recessions  at  the 
higher  elevations  are  more  reliable  indices  of  shore  retreat. 


Distance  from  Baseline  (meters) 


Fig.  7 Adjust- 
ments of  the 
Upper  Profile 
at  Station  5. 


Average  recession  (1967  to  1971)  of  the  176.30  m shoreline  totaled 
6.5  m.  For  this  four  year  period  transgression  was  responsible  for  more 
than  50%  of  the  total  shoreline  retreat. 

b.  Nearshore.  The  nearshore  profile  is  dominated,  from  near  water's 
edge  to  a point  approximately  500  m from  shore,  by  a sequence  of  from  four 
to  five  longshore  bars.  These  bars  are  persistent  year-round  features,  but 
are  not  stationary  (Figure  8).  On  the  north  side  of  the  harbor  where  four 
bars  were  persistent  from  year  to  year  throughout  the  study,  the  inner 
three  bars  migrated  an  average  of  26  m toward  the  shore,  and  rose  in  eleva- 
tion about  0.5  m during  the  same  four  year  period.  The  cross-sectional 
geometries,  aerial  relationships,  and  migration  of  multiple  longshore  bars 
are  discussed  in  detail  elsewhere  for  a larger  region  encompassing  Pent- 
water  Harbor  (Hands,  1976). 

4.  Interpretation. 

The  120  year  rate  of  historic  recession  for  a typical  stretch  of  uncon- 
solidated lake  shore  is  about  0.37  m/yr  (1.2  ft/yr) . Rates  of  recession  are 
not  however  constant;  periods  of  accelerated  recession  occur  during  years 
of  high  lake  level.  If  measurements  of  recession  obtained  during  the  recent 
episode  of  high  water  are  divided  into  two  nearly  equal  time  intervals 
(1967  to  1969,  1969  to  1971),  each  reflecting  equal  submergence  (0.2  m) , 
then  recession  of  the  highest  common  ihofiziAne  (176.30  m)  was  about  the  same 
for  both  periods  and  totaled  6.5  m. 

Total  shoreline  retreat  exceeded  recession  by  a factor  of  more  than 
two.  The  difference  between  the  total  retreat  (14.6  m)  and  recession  (6.5) 
is  transgression  (8.1).  In  other  words,  in  addition  to  6.5  m lost  by  ero- 
sion, 8.1  m of  shore  has  been  lost  by  submergence  beneath  the  elevated  lake 
levels . 

Figure  9 shows  the  total  changes  in  average  position  of  bar  crests,  bar 
troughs,  and  the  176.3  ihofi££Anz  between  1967  and  1971.  Changes  in  eleva- 
tion of  crests,  troughs,  and  ikoAeXtne  were  essentially  the  same  (0.55, 

0.47,  and  0.51  m respectively).  Average  horizontal  changes  were  25  m for 
the  crests,  24  m for  the  troughs,  but  only  6.5  m for  the  1967  ihomtAm. 

An  important  question  is  whether  the  documented  shoreline  recession 
represents  full  or  only  partial  adjustment  to  submergence.  Per  Bruun  (1962) 
has  hypothesized  that  there  is  an  equilibrium  form  which  beaches  tend  to 
adher  to,  and  if  sea-level  rises,  the  equilibrium  form  will  be  shifted  up- 
ward and  landward.  Lake  Michigan  bars  moved  up  by  an  amount  equal  to  the 
mean  rise  in  lake  level  during  the  four  year  period.  However,  while  main- 
taining a fixed  depth,  the  bars  encroached  on  the  shoreline  (Figures  8 and  9) . 
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Figure  8.  Landward  Migration  of 
Longshore  Bars  Over  a Four 
Year  Period  of  Rising  Lake 
Levels . 


If  the  1967  profiles  more  nearly  approximated  the  equilibrium  form  than  did 
the  later,  steeper  profiles,  then  considerable  additional  retreat  must 
occur  after  lake  levels  stop  rising,  to  flatten  the  profiles  back  to  their 
1967  configuration.  The  difference  in  bar  and  shoreline  migration  (25-15  = 
10  m)  is  interpreted  as  a lag  in  the  recession  of  the  shoreface  or  upper 
beach.  In  order  to  adjust  completely  to  the  elevated  lake  surface,  the 
upper  part  of  the  profile  will  probably  recede  farther  landward  by  continued 
rapid  erosion  until  it  has  increased,  by  roughly  10  m,  its  separation  from 
the  nearshore  bars.  Assuming  shore  erosion  supplies  a volume  of  sediment 
sufficient  to  readjust  the  nearshore  profile,  a crude  sediment  budget  was 
calculated  (Hands,  1975) . The  results  were  in  substantial  agreement  with 
the  preceding  prediction  inasmuch  as  mass  balance  provided  an  even  larger 
estimate  of  the  additional  retreat  necessary  to  re-establish  equilibrium; 
the  final  ratio  of  recession  to  submergence  was  estimated  to  be  about  60:1. 


Figure  9.  Migration  of  Bars  and  Shoreline  (1967  to  1971).  Migration  of  Bar 
Crests,  D , Trough  Thalwegs,  t , and  the  Water's  Edge,  o , from  their  Mean 
Positions  in  1967  to  Their  Mean  Positions  in  1971.  Based  on  Profile 
Measurements  at  the  Three  Stations  North  of  Pentwater  Harbor. 


Recently  collected  profile  data  provide  more  extensive  coverage  both  along 
shore  and  offshore,  and  will  thus  provide  a better  basis  for  future  refine- 
ments of  the  sediment  budget  . 


Ill  USE  OF  DIFFERENTIAL  UPLIFT  TO  MODEL  SUBSIDENCE 

1 . Regional  Trends  In  Recession  Rates. 

By  selecting  from  Powers'  report  those  stations  initially  surveyed 
between  1830  and  1838  (94  in  number)  and  then  averaging  them  by  county, 
evidence  for  an  unreported  regional  trend  in  long  term  recession  rates  is 
obtained  (Fig.  10).  The  relatively  large  variation  in  recession  rate 
encountered  as  one  moves  from  one  station  to  the  next  along  the  shore, 
cautions  against  putting  too  much  confidence  in  the  mean  rate  derived  for 
any  given  county.  An  average  of  only  four  sample  rates  isn't  a very  stable 
estimate  of  the  true  mean  rate,  and  the  number  of  measurements  per  county 
varies.  Yet  recurrent  increases  in  recession  rates  toward  the  south  on 
both  sides  of  the  lake  support  the  hypothesis  that  a regional  trend  actual- 
ly exists.  Evidence  of  this  trend  has  not  been  found  in  any  other  data  on 
Lake  Michigan  shore  recession,  but  no  other  data  set  has  an  area  coverage 
and  time  span  comparable  to  Powers'.  The  fact  that  Powers  did  not  suggest 
a regional  trend,  and  was  apparently  unaware  of  the  evidence  for  one  in  his 
data,  eliminates  the  possibility  of  even  unconscious  bias  in  the  reduction 
of  survey  notes  and  compilation  of  recession  rates.  Hhat  the.i  could  be 
responsible  for  this  regional  trend  in  recession  rates? 

2 . Possible  Explanations. 

The  rate  of  bluff  recession  depends  basically  on:  a)  lake  level  be- 
havior, b)  erosive  forces,  c)  resistance  of  the  shore  deposit,  and  d)  the 
offshore  profile.  Waves  are  the  primary  source  of  energy  needed  to  do  the 
work  of  shore  erosion,  so  the  possibility  that  fetch  and  wind  conditions 
might  give  rise  to  a trend  in  wave  energy  was  examined.  Visual  wave  obser- 
vations have  been  taken  for  a number  of  years  at  various  stations  around 
Lake  Michigan  (LEO  Data,  unpubl.).  The  means  of  observed  breaker  heights 
for  a three  year  period  common  to  24  stations  were  used  to  estimate  the 
average  wave  power  entering  the  surf  zone  at  each  station.  Assuming  the 
wave  power  delivered  to  the  surf  zone  is  proportional  to  the  5/2  power  of 
the  breaker  height  (CERC,  1973)  , 
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and  furthermore  assuming  that  the  power  of  the  average  observed  breaker  is 
a good  index  of  average  breaker  power,  the  distribution  of  Hj,  5/2  around  the 
lake  was  plotted  (Fig.  11).  Neither  raw  Hj,,  Hb  nor  an  index  of  breaker 

power  gradient,  AHb  (where  AX  is  the  shoreline  distance  between 

stations)  suggested  any  regional  trend  in  erosive  forces. 

The  possibility  that  northern  shores  might  benefit  from  a longer  period 
of  isolation  from  winter  storm  waves  by  pack  ice,  was  eliminated  by  consult- 
ing an  ice  atlas  (Rondy,  1969). 

Essentially  all  the  eastern  shore  of  the  lake  can  be  described  as  alter- 
nating sections  of  sand  dunes  intersected  by  glacial  moraines.  Nothing  in 
the  distribution  of  shore  types,  height  of  the  sand  dunes,  or  morainal  bluffs 
(Hands,  1970)  suggest  any  regional  trend  in  the  shores'  resistance  to  erosion. 
Although  nearshore  slopes  do  gradually  flatten  at  the  southern  end  of  the 
lake,  if  the  nearshore  slope  had  any  effect  on  recession  rates,  it  would 
tend  to  decrease  recession  in  the  southern  area  rather  than  the  reverse, 
which  was  observed. 

The  last  independent  variable  considered  is  lake  level.  Annual  and 
long  term  fluctuations  in  lake  level  affect  all  parts  of  the  basin  equally. 


Figure  10.  Rates  of  bluff  reces- 
sion [m/yr]  averaged  by  counties. 
Tic  marks  indicate  number  and  lo- 
cation of  measurements.  *For 
Porter  Co.,  the  average  is  for 
the  period  1927  to  1957.  All 
other  averages  are  for  the 
period  1830's  to  1950's. 


Figure  11.  Distribution  of 
"average  wave  power"  based  on 
visual  observations  of  breaker 
height  at  24  LEO*  stations 
each  reporting  three  years  of 
data.  *Littoral  Environment 
Observation  Program  at  the 
U.S.  Army,  Coastal  Engineering 
Research  Center. 


Long  term  relative  submergence  would,  thus,  have  to  be  a function  of 
differential  crustal  motion. 

3.  Evidence  of  Differential  Uplift  of  the  Lake  Michigan  Basin. 

Abandoned  strandlines,  relict  from  ancestral  lakes,  are  particularly 
well  developed  in  the  Great  Lakes  basin  above  the  elevation  of  present  lake 
surfaces.  Early  geologists  traced  these  abandoned  shore  features  for 
hundreds  of  kilometers  and  found  that  they  were  not  level,  but  rose  in  ele- 
vation toward  the  north.  Assuming  a given  strandline  formed  approximately 
synchronously  along  its  length,  and  was  initially  close  to  level,  the  ob- 
served tilt  was  Interpreted  as  a measure  of  crustal  uplift  (Fig.  12) . For 
Instance  in  the  northern  Lake  Michigan  basin  the  Niplssing  shoreline  rises 
about  12  m in  200  km,  indicating  12  m of  differential  uplift  occurred  since 
these  features  originally  formed. 

Water  level  records  form  a second  and  independent  source  of  evidence 
for  crustal  motion.  At  a number  of  harbors  on  Lake  Michigan,  water  level 
records  extend  back  past  the  turn  of  the  century.  Differences  in  lake  sur- 
face elevations  recorded  at  separate  locations  vary  with  time,  but  the  time 
series  of  differences  should  have  a stationary  mean  if  the  two  locations  are 
not  subject  to  differential  crustal  movement.  Analysis  of  crustal  movement 
using  this  approach  goes  back  to  Gilbert  (1898) . The  total  record  of  lake 
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Figure  12.  Geologic  Evidence  of  Basin  Tilt.  Elevation  of  Ancient  (3500  yr. 
B.P.)  Strandline  Features  in  Feet  from  Leverett  and  Taylor  (1915). 


Figure  13.  Regression  of  Lake  Level  Differences  at  Different  Gages  Indicate 
Basin  Tilt  (m/100  yrs);  from  Kite  (1972). 
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level  data  has  been  recently  reviewed,  edited,  adjusted,  and  analyzed  by 
Kite  (1972).  His  estimate  of  basin  tilt,  shown  in  Figure  13,  is  based  on 
lake  level  records  ranging  from  35  to  110  years  in  length.  Differences  in 
water  levels  recorded  at  Milwaukee  and  Sturgeon  Bay  for  example.  Increase 
over  the  65  year  period  of  common  record.  The  Increase  is  approximately 
linear  with  time,  and  at  a rate  of  2 mm/yr. 

A third  Independent  line  of  evidence  for  crustal  movement  comes  from 
geodetic  leveling.  By  comparing  adjusted  data  obtained  in  1929  and  1955 
level  surveys,  Holdahl  produced  the  unpublished  map  shown  in  Figure  14. 
Differences  in  measured  elevations  were  found  to  Increase  along  the  level 
line  running  northwesterly  from  Chicago.  By  extrapolation  the  Lake 
Michigan  basin  is  thought  to  have  tilted  0.18  m in  the  26  years  between 
surveys.  Between  Milwaukee  and  Sturgeon  Bay  this  would  again  be  roughly 
equivalent  to  2 mm/yr  of  differential  uplift.  These  results  were  described 

CRUSTAL  MOVEMENT  BASED  ON 


Figure  14.  Comparison  of  First  Order  Level  Net  of  1929  with  First  Order 
Releveling  in  1955  (Meade,  1972)  indicates  Basin  Tilt. 

(Meade,  1971)  as  tentative  and  subject  to  revision,  but  for  present  purposes 
they  are  more  than  adequate;  the  general  pattern  of  an  active  regional  tilt 
in  the  Lake  Michigan  basin  is  confirmed  by  the  similarity  of  results  from 
independent  lines  of  Investigation. 

4 . Effect  of  Uplift  on  Long  Term  Recession  Rates. 

Could  the  relative  submergence  of  the  southern  end  of  the  Lake  Michigan 
basin  explain  the  greater  bluff  recession  that  has  occurred  there?  It  is 
usually  assumed  that  even  if  crustal  uplift  is  an  active  process,  it  is  too 
slow  to  be  a significant  factor  in  contemporary  erosion  problems.  Given  the 
uncertainties  in  computed  rates  of  uplift  and  bluff  recession,  a possible 
relationship  between  the  two  should  be  examined  in  the  simplest  manner  pos- 
sible. In  order  to  compare  them  quantitatively,  both  will  be  approximated 
by  their  linear  trends.  This  is  not  meant  to  imply  that  recession  is  stric- 
tly a function  of  lakeshore  position.  As  pointed  out  earlier  there  are  many 
factors  affecting  recession.  The  attempt  here  is  merely  to  obtain  a quanti- 
tative estimate  of  how  subsidence,  taken  by  itself,  affects  recession.  The 
regression  coefficient  is  a good  estimator  of  this  effect  as  none  of  the 
other  factors  is  thought  to  exert  a regional  control.  The  remaining  scatter 
in  Figure  15  illustrates  the  combined  effect  of  these  other  variables,  which 
as  expected  is  considerable.  The  least  square  regression  coefficient  is  19 
+10  (x  + 2s)  m per  century  per  100  km  along  the  lake  axis;  this  trend  is 


statistically  significant  even  at  the  1%  level.  From  Figure  15  It  seems 
that  this  trend  may  be  an  expression  of  a northerly  decrease  In  an  upper 
bound  on  recession,  with  low  recession  values  distributed  fairly  uniformly 
up  and  down  the  lakeshore.  IVo  low  values  near  60  and  80  km  contribute 
heavily  to  this  Impression.  It  Is  speculated  that  greater  Investment  In 
shore  protection  may  explain  the  occurrence  of  some  low  recession  values  In 
the  earlier  settled  southern  portion  of  the  basin. 

The  period  of  time  covered  by  the  water  level  records  and  by  the  bluff 
recession  data  discussed  In  Section  1 are  roughly  the  same.  The  tilt  of 
the  basin  estimated  from  Figure  13  and  14  would  be  .063  and  .087  m per  cen- 
tury per  100  km  along  the  axis  of  the  lake.  Thus  each  centimeter  of  subsi- 
dence caused  somewhere  between  one  and  four  meters  of  recession  If  the  trend 
In  recession  Is  to  be  attributed  solely  to  subsidence. 

It  Is  noted  that  the  longshore  trend  In  recession  rates  Is  not  so  pro- 
nounced on  the  east  shore  as  on  the  west.  Nevertheless  the  hypothesis  that 
the  apparent  trend  might  have  arisen  by  chance  (H^  :S  - 0)  would  still  be 
rejected;  under  standard  assumptions  even  If  the  test  were  applied  to  the 
east  coast  data  taken  by  Itself,  t = 2.4,  n = 33).  With  respect  to  the 
smaller  Increase  In  recession  along  the  east  shore.  It  Is  Interesting  to 
note  that  the  east  shore  Is  backed  by  high  dunes  and  glacial  deposits  rich 
In  sand  and  gravel.  Each  meter  of  bluff  retreat  on  the  east  shore  would 
supply  a greater  volume  of  beach  material  than  on  the  typically  lower 
western  shore . 


Figure  15.  Long  Term  Rate  of  Recession  Versus  Position  Measured  Along 
Central  Axis  of  the  Lake.  Values  from  the  West  Shore  are  Shown  as  • , 
East  Shore  as  t . 
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IV 


CONCLUDING  REMARKS 


Subsidence  Increases  the  rate  of  erosion  In  unconsolidated  deposits  by 
allowing  waves  to  reach  bluffs  that  were  formerly  above  the  elevation  of 
wave  uprush.  Moreover,  waves  lose  less  energy  In  passing  over  submerged 
offshore  shoals.  The  Increased  erosion  can  be  viewed  as  an  adjustment  of 
the  beach  to  new  conditions  imposed  by  subsidence.  The  quantitative  re- 
lationship between  subsidence  and  the  rate  of  shore  retreat  has  however, 
received  little  attention.  From  the  standpoint  of  shore  protection  there 
are  many  other  over-riding  variables  that  effect  erosion.  And  In  areas  of 
extreme  subsidence,  there  are  many  other  damaging  consequences:  water 
supply,  structural  failures,  etc.  This  study  demonstrates  that  subsidence 
can  have  a measurable  effect  on  shore  retreat.  The  effect  could  have  pro- 
found local  impact  - consider  the  resulting  decrease  in  property  value  if 
a formerly  stable  barrier  Island  where  subjected  to  moderate  subsidence 
over  a fifty  year  period. 

The  shore  of  Lake  Michigan  retreated  15  m during  a four  year  period  in 
response  to  a 0.5  m Increase  in  mean  water  elevation,  and  because  a broad 
section  of  the  submerged  profile  responded  by  moving  25  m landward,  it 
appears  that  the  shore  was  lagging  behind  in  its  adjustment  to  high  lake 
levels.  Total  recession  required  to  reestablish  equilibrium  after  a 0.5  m 
coastal  submergence  in  four  years,  extends  years  beyond  the  period  of  sub- 
sidence. Based  on  documented  lake  shore  retreat,  profile  change,  and  a 
rough  balancing  of  the  sediment  budget,  it  is  estimated  that  25  to  30  m of 
recession  would  be  required  to  readjust  the  shore  to  the  0.5  m rise  in  lake 
level.  This  gives  a ratio  of  expected  recession  to  submergence  of  approxi- 
mately 60  to  1. 

The  amount  of  shore  retreat  that  would  result  from  a much  slower  rate 
of  subsidence  was  interpreted  using  records  of  long  term  bluff  retreat  and 
crustal  motion.  The  Lake  Michigan  basin  has  been  tilting  upward  toward  the 
north.  The  120-year  mean  recession  rates  show  a similar  trend;  with  reces- 
sion rates  increasing  in  the  direction  of  greatest  subsidence.  Regional 
tilting  seems  to  be  the  best  explanation  for  the  regional  variation  of 
recession.  Under  this  assumption  each  cm  of  slow  submergence  would  be 
responsible  for  from  one  to  four  m of  shore  recession;  the  ratio  of  recession 
to  submergence  would  be  between  100:1  and  400:1. 

Shore  retreat  is  a non-linear,  time -dependent  function,  and  so  apparen- 
tly is  the  relationship  between  recession  and  subsidence  rates.  Great  Lakes 
studies  suggest,  that  long  term  subsidence  may  cause  several  times  the  reces- 
sion that  would  result  during  a short  period  of  equal,  but  rapid  subsidence. 

It  takes  a number  of  years  for  the  beach  profiles  to  equilibrate.  Moreover, 
slow,  long  term  water  level  adjustments  may  permit  littoral  forces  to  spread 
shore  eroded  material  across  a wider  submerged  profile;  and  therefore,  in  the 
long  term  require  a greater  volume  of  shore  eroded  material  to  adjust  to  sub- 
me  rgen  ce . 

Coastal  areas  with  geomorphology,  geology,  and  wave  exposure  similar  to 
the  study  area  may  be  expected  to  recede  at  rates  roughly  on  the  order  of 
those  measured  if  subjected  to  the  same  conditions  of  subsidence.  Wherever 
possible,  shore  recession  caused  by  subsidence  should  be  determined  along 
with  measurements  of  subsidence.  If  enough  additional  data  can  be  acquired, 
it  may  be  feasible  to  establish  some  functional  relationship  between  subsi- 
dence and  shore  recession  that  will  be  valid  for  a range  of  subsidence  rates. 
Resulting  relationships  may  then  apply  to  broad  classes  of  coastal  conditions. 
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